Interaction-induced delocalization of quasiparticle pairs in the Anderson insulator 
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It has recently been established that a short-range interaction can strongly delocalize a pair 
of particles moving in a disordered potential. We investigate whether an analogous effect 
exists also for pairs of quasiparticles in the Anderson insulator by employing an approximate 
numerical evaluation of the two-particle Green function of the many-body system. At the 
Fermi energy, the quasiparticle pair is localized on the scale of the single-particle localiza- 
tion length. The delocalization effect is recovered when the excitation energy of the pair is 
comparable to the bandwidth. 

PACS numbers: 72.15.R, 71.30. 



The localization of particles by a random potential 
has been studied extensively over the past decades [Q. 
While wavefunctions are exponentially localized in one 
and two dimensions even for arbitrarily weak disorder, 
there exists a transition between extended and localized 
states in three dimensions. A non-interacting degenerate 
Fermi gas in a disordered potential is conducting when 
the Fermi energy lies in the region of extended states 
(metallic phase) while it is insulating when the Fermi 
energy lies in the region of localized states (Anderson 
insulator). The situation becomes much more complex 
and controversial once interactions between the particles 
are included. Interactions may lead to a metal-insulator 
transition (Mott transition) even in the absence of dis- 
order; they also result in rich magnetic behavior j^]. In- 
terest in the interplay of disorder and interactions was 
recently renewed in part by experiments measuring per- 
sistent currents in mesoscopic normal-metal rings whose 
amplitude could not be explained by theories for non- 
interacting electrons j|] . 

An original approach to studying localization in the 
presence of interactions was recently taken by Shepelyan- 
sky 01 and previously by Dorokhov ||] . They considered 
the much simpler problem of two interacting particles in 
a random potential and predicted that the interaction 
can lead to a significant delocalization of the pair. She- 
pelyansky Q and Imry Q speculated that an analogous 
effect could exist for pairs of quasiparticles at finite par- 
ticle density which may have significant consequences for 
the thermodynamic and transport properties of the An- 
derson insulator. It is the purpose of the present paper to 
study the propagation of quasiparticle pairs in the Ander- 
son insulator in the presence of a short-range interaction 
within a numerically tractable approximation to the two- 
particle Green function. Our approximation is motivated 
by both physical and analytical considerations. 

We briefly recall some of the pertinent results for 
the two-particle problem. Shepelyansky studied two 
particles in a random potential interacting by a short- 
range interaction. Whenever the two particles are local- 



ized far apart compared to the single-particle localization 
length ^1, the effect of the interaction is only exponen- 
tially small. By contrast, the two particles can prop- 
agate as a pair over a distance ^2 which can be much 
larger than when they are localized within about one 
single-particle localization length of each other. Specif- 
ically, Shepelyansky predicted by an approximate map- 
ping to a banded-random-matrix model which he stud- 
ied numerically that the two-particle localization length 
^2 satisfies ^2/^1 — (Ci/32)(u/t)^, where u denotes the 
disorder strength and t the bandwidth (hopping matrix 
element). The effect is essentially independent of the 
sign of the interaction and the statistics of the parti- 
cles. Subsequently, the existence of this effect was con- 
firmed using both a Thouless-type scaling argument j|] 
and numerical calculations . Numerical work ^ com- 
puting ^2 directly from the Green function of the two- 
particle problem showed that in strictly one dimension 
the two-particle localization length at the center of the 
band satisfies a scaling relation ^2/Ci = /(^Ci/i) with 
f{x) ~ 1/2-1- C\x\ which, while qualitatively confirm- 
ing the effect, is inconsistent with the original prediction 
mentioned above (C is a numerical constant). 

The most important additional feature of the degen- 
erate Fermi system is the Pauli principle. Two parti- 
cles propagating in a disordered potential sufficiently far 
from the band edges have ample phase space for scat- 
tering. On the other hand, phase space for scattering is 
severely restricted at finite density since all states below 
the Fermi energy are effectively blocked. One may ex- 
pect that this leads to a suppression of the delocalization 
at low pair excitation energy. With increasing excitation 
energy, the phase-space restriction becomes less severe 
and the delocalization effect should be recovered. We 
study the influence of the Pauli principle on the prop- 
agation of quasiparticle pairs within an approximation 
analogous to the Cooper problem of superconductivity: 
We consider two particles above the Fermi energy which 
interact with each other, but not with the electrons in the 
Fermi sea except via the exclusion principle. The Fermi 
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FIG. 1. Examples of elementary second-order diagrams 
in the particle-particle channel with time running upwards. 
While diagram (a) contributes to leading order in the 
low-density limit, diagrams (b) and (c) are suppressed since 
they involve particle-hole-pair creation and annihilation. 

sea is defined by filling up single-particle eigenstates of 
the disorder potential below the Fermi energy Ep. In or- 
der to avoid the superconducting instability, we consider 
only repulsive interactions. 

One can motivate this approximation also by a quite 
different line of thought. It is plausible that the two- 
particle effect is most likely to survive at low particle 
density. When considering a short-range interaction in 
three dimensions, only ladder diagrams involving no in- 
termediate hole excitations contribute to leading order 
in the low-density limit Hence, the Cooper prob- 
lem is in fact equivalent to a low-density approximation. 
This can be understood physically by considering some 
second-order diagrams as shown in Fig. 1. While the 
ladder diagram in Fig. 1(a) involves only particle exci- 
tations at the intermediate stage, the diagrams 1(b) and 
1(c) involve intermediate hole excitations. The latter di- 
agrams are suppressed at low density because they in- 
volve particle- hole-pair creation. For these reasons, we 
will only study particle-particle excitations in this paper. 
We expect that our results apply qualitatively also to 
hole-hole and particle- hole excitations. 

We now proceed to work out the Cooper approxima- 
tion within the diagrammatic approach. We consider 
spinless fermions on a one-dimensional lattice with re- 
pulsive nearest-neighbor interaction of strength u and 
subject to a random potential, 

n n 

+'"^atal+ia„+ia„ . (1) 
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Here, a„ denotes the fermionic annihilation operator of 
a particle at site n, and w„ G [—W/2,W/2] denotes the 
random on-site energies. In the following, energies will 
be measured in units of the hopping matrix element t and 
lengths in units of the lattice spacing a. As motivated 
above, we focus on the retarded two-particle Green func- 
tion in the particle-particle channel. Analogous to our 
approach j|] to the two-particle problem, we consider 
only processes where the quasiparticles are created on 



neighboring sites (m, m+1) at time t = and destroyed 
at time r on a different set of neighboring sites (n, n + 1). 
The corresponding amplitude is given by 

FPP(n,m;T) = -z0(T)(O|a„(T)a„+i(T)at„+i(O)a1„(O)|O) 

— FPP{72,m;E)exp{-iET/h). (2) 
-oo 27r 

The localization length for coherent propagation of 
the quasiparticle pair can be defined by the exponential 
decrease of these matrix elements with distance, 
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where (...) denotes a disorder average. When consider- 
ing the nearest-neighbor interaction in (^), a closed equa- 
tion can be derived for the matrix elements i^PP(n, to; E) 
within the ladder approximation. It is this enormous sim- 
plification which makes our numerical calculations feasi- 
ble. 

Computing _FPP(n, to; E) for non- interacting fermions, 
one readily obtains 



F^''{n,m;E) 



(4) 



Here, denotes the single-particle energy of the eigen- 
state (t>i{n) of the disorder potential, i] is a positive in- 
finitesimal, and Bij{n) = (j)i{n)(l)j{n+l) — 4>i{n+l)4>j{n). 
Summing the ladder diagrams and considering only the 
contributions of intermediate particle-particle states, one 
obtains the two-particle Green function in the Cooper ap- 
proximation, 
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where Fpp{E) is viewed as a matrix in site space. It is 
instructive to compare these equations to those for the 
two-particle problem [|j . The difference between the two 
cases is that the sum in (^ is only over states above 
the Fermi energy while the corresponding equation for 
the two-particle case involved an unrestricted summa- 
tion over all states. It is possible to extend our ap- 
proach to the full ladder approximation. Then, Eq. (^) 
would also contain a summation over hole-hole excita- 
tions. However, we decided to use the Cooper approx- 
imation because it can be motivated systematically in 
the low-density limit. We do not expect that our re- 
sults change qualitatively when hole-hole excitations are 
included. 

The first factor on the right-hand side of (|^) decreases 
on the scale of the single-particle localization length. 
Hence, we focus on the second factor in the following 
from which any long-range behavior must arise. For nu- 
merical purposes, it is instructive to interpret this term 
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as the Green function for the "Hamiltonian" Fq^{E) at 
"energy" 1/u. Since Fq^{E) is a banded matrix whose 
bandwidth is of the order of the single-particle local- 
ization length, this allows us to compute the localiza- 
tion length S,q by employing the efficient recursive Green- 
function technique for banded Hamiltonian matrices . 

For our numerical calculations we chose a Fermi energy 
Ep = — 1, corresponding roughly to quarter filling. We 
have studied the localization length of quasiparticle pairs 
for nine values of the pair energy E ranging from E = —2 
to -E = in steps of AE = 0.25. This corresponds to 
excitation energies e = E ~ 2Ep of the quasiparticle pair 
between e = and e = 2. For each value of the excitation 
energy we used three values of disorder W = 1.5,2, and 
3, corresponding to single-particle localization lengths at 
the band center — 46.7,26.2, and 11.7, respectively, 
and five values of the interaction u ~ 0.1,0.3,0.5,0.75 
and 1.0. For comparison, we have also computed ^2 of 
the two-fermion problem for the same set of parameters. 
The computations were done for lattices with 500 sites, 
and we averaged over 50 realizations of the disorder. We 
briefly comment on our choice of Ep. As argued above, 
our approximation can be motivated by the low-density 
limit. In view of this, our choice oi Ep may seem quite 
large. We chose this somewhat larger value of Ep since 
the two-particle effect becomes increasingly weaker as one 
moves away from the center of the band, mostly because 
the single-particle localization length decreases. Hence, 
for smaller values of Ep, it would be difficult to identify 
a suppression of the delocalization effect relative to the 
two-particle effect. 

In Fig. 2 we exhibit our results for the two pair energies 
E = —2 and E ~ corresponding to excitation energies 
e = and e = 2, respectively. We plot ^g{E)/Ci{E/2) 
as a function of u£_i{E/2)/t aiialogous to the scaling plot 
for the two-particle problem.cl While we do not find the 
same scaling any more in the many-particle problem, this 
still allows us to exhibit data for all values of the interac- 
tion and the disorder strength in a single plot. Clearly, 
the delocalization effect, disappears almost completely for 
zero excitation energy.!!! On the other hand, the delocal- 
ization of the quasiparticle pair is close to that of the 
two-particle problem for e = 2 as seen by comparison 
with the two-particle result at £^ = also plotted in 
Fig. 2. 



*Note that the single-particle localization energy should be 
evaluated at half the pair energy. We generated the single- 
particle localization lengths for these energies numerically by 
the recursive Green function method [|_0| using systems with 
10® sites. 

^Note that we find ^2 < Ci/2 for e = 0. This is a consequence 
of neglecting the first factor in Eq. However, there is 

indeed no delocalization effect in this case, because the range 
of both factors in (H) is independent of u and W . 
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FIG. 2. Plot of ^g{E)/^i{E/2) vs <i(_B/2)/f for two dif- 
ferent excitation energies. Five values of u are included for 
each of the three values of disorder W = 3 (triangles) , W = 2 
(crosses), and W = 1.5 (diamonds). The two-fermion effect 
is also shown for comparison (dashed line). Evidently, the 
delocalization effect has disappeared at the Fermi energy and 
is recovered with increasing excitation energy. One also ob- 
serves deviations from the scaling behavior of the two-fermion 
case. 

In view of the absence of the delocalization effect at 
zero-excitation energy, it is important to identify the rel- 
evant scale on which the effect is recovered with increas- 
ing excitation energy. To investigate this question, we 
have considered the enhancement of relative to the 
enhancement of ^2 for the two-particle problem as mea- 
sured by 
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e2(ti)-6(w = o.i) ■ 

We have computed this quantity as a function of the ex- 
citation energy e for different disorder and interaction 
strengths. Between e = and e = 2, the quantity should 
vary roughly between zero and one. Figs. 3(a) and (b) 
show that Ru,w[^ at fixed u is independent of the dis- 
order strength VF.cl Moreover, a comparison of Figs. 3(a) 
and (b) indicates that Ru,w{f-) is also independent of the 
interaction strength u. This implies that the bandwidth 
t is the relevant scale on which the delocalization effect 
is recovered with increasing excitation energy. This re- 
sult is in qualitative agreement with the prediction from 
Thouless-type arguments by Imry 

Our results were obtained using a simple approxima- 
tion to the two-particle Green function. Even within im- 



*The accuracy of the data for -Ru.w.(e) is limited because 
both numerator and denominator in (H) are small numbers. 
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FIG. 3. Plot of the delocalization effect Ru,w{^) as defined 
in Eq. (^) as a function of excitation energy for two values 
of u and three values of disorder W — 3 (triangles) , W — 2 
(crosses), and W = 1.5 (diamonds). Note that Ru,w{e) is 
insensitive to both disorder and interaction strength. 
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proved approximations we do not expect a delocaliza- 
tion effect analogous to that for two particles to reap- 
pear at low excitation energies because the suppression 
of the two-particle effect at finite density is entirely due 
to the Pauli principle. This expectation is also in accord 
with standard phase-space arguments for the degenerate 
Fermi system. On the other hand, our approximation 
does not allow us to exclude the possibility that the local- 
ization properties of the Anderson insulator are strongly 
affected by true correlation effects. 

We conclude by discussing the physical implications of 
our results. For thermodynamic and low- frequency trans- 
port properties in linear response, the relevant excitation 
energies are of the order of the temperature. In these 
cases, the maximal delocalization of quasiparticle pairs 
is not expected to be observable since the bandwidth 
will generally be much larger than the temperatures at 
which quantum-interference effects such as Anderson lo- 
calization are observable. In this sense, the quasiparti- 
cle delocalization at finite density is significantly weaker 
than the two-particle effect. However, it may still be pos- 
sible to observe the incipient delocalization at excitation 
energies which are low compared to the bandwidth. In 
particular, this may be possible in higher dimensions ||] 
where the delocalization effect is expected to be more pro- 
nounced and where the single-particle localization length 
can become exceedingly large (e.g., close to the Ander- 
son transition). Finally, we note that there may also be a 
possibility of detecting the delocalization of quasiparticle 
pairs in high-frequency experiments. 

We enjoyed helpful discussions with G. Barkema, J. 
Miiller, and H.A. Weidenmiiller. 
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